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Abstract 

We have searched for the decay K$ — > 37r° with the KLOE experiment at DA<I>NE 
using data from e + e~ collisions at a center of mass energy W ~ for an inte- 
grated luminosity L = 450 pb _1 . The search has been performed with a pure K$ 
beam obtained by tagging with Kl interactions in the calorimeter and detecting six 
photons. We find an upper limit for the branching ratio of 1.2 x 10~ 7 at 90% C.L. 

Key words: e+e' collisions, DA$NE, KLOE, rare K s decays, CP, CPT 



1 Introduction 

The decay K s — > 3n° violates CP invariance. The parameter r] 000 , defined as 
the ratio of K s to K L decay amplitudes, can be written as: r] 000 = A(K S — > 
3tt°)/A(Kl — > 37T°) = e + eg 00 , where e quantifies the K$ CP impurity and 6q 00 
is due to a direct CP- violating term. Since we expect eQ 00 e [1], it follows 
that r?ooo~e. In the Standard Model, therefore, ~BR(K S -> 3tt°) ~1.9 x 10~ 9 
to an accuracy of a few %, making the direct observation of this decay quite 
a challenge. 

The best upper limit on BR(K S — > 3n°) from a search for the decay 
was obtained by the SND experiment at Novosibirsk. They find BK(K,s — > 
3tt ) < 1.4 x 10~ 5 at 90% C.L. [2]. CPLEAR has pioneered the method of 
searching for interference between K$ and Kl decays. Interference results 
in the appearance of a term 3? (?7ooo) cos(Amt) — 3 (^ooo) sin(Amt) in the 
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decay intensity. K (7/000) an d S (77000) are obtained from a fit, without dis- 
criminating between Kl or K$ — > 37r° decays. In this way CPLEAR finds 
r/ 000 = (0.18±0.15)+i (0.15±0.20) [3]. The NA48 collaboration [4] has recently 
reached much higher sensitivity. By fitting the K s jK L — > 3ir° interference pat- 
tern at small decay times, they find K (77000) = —0.002 ± 0.011 sta t ± 0.015 sys 
and 3= (77000) = —0.003 ± 0.013 sta t ± 0.017 sys , corresponding to BR(K S — > 
3tt°) < 7.4 x 1(T 7 at 90% C.L. The sensitivity to CPT violation via unitar- 
ity [5] is now limited by the error in 7/ H = A{Kl — > tt + tt~)/A(Ks — > 7r + 7r~). 

We report in the following an improved limit from a direct search for the 
37T° decays of the K$- Apart from the interest in confirming the Standard 
Model, knowledge of 7/000 allows tests of the validity of CPT invariance using 
unitarity. 



2 DA$NE and KLOE 



The data were collected with the KLOE detector [6-9] at DA$NE [10], the 
Frascati factory. DA$NE is an e + e~ collider operated at a center-of-mass 
energy W ~ 1020 MeV, the mass of the meson. Positron and electron beams 
of equal energy collide at an angle of tt— 0.025 rad, producing mesons nearly 
at rest (p<p ~ 12.5 MeV). mesons decay 34% of the time into nearly collinear 
pairs. Because J PC (0) = 1 , the kaon pair is in a C-odd antisymmetric 
state, so that the final state is always Ks-Kl- Detection of a Kl signals the 
presence of a Ks of known momentum and direction. We say that detection 
of a K L "tags" the K s . 

The KLOE detector consists of a large cylindrical drift chamber (DC), 
surrounded by a lead/scintillating- fiber electromagnetic calorimeter (EMC). 
A superconducting coil around the calorimeter provides a 0.52 T field. The 
drift chamber, 4 m in diameter and 3.3 m long, is described in Ref. 6. The mo- 
mentum resolution is o"(jo_i_)/p_i_ ~ 0.4%. Two track vertices are reconstructed 
with a spatial resolution of ~ 3 mm. The calorimeter, described in Ref. 7, is di- 
vided into a barrel and two endcaps, for a total of 88 modules, and covers 98% 
of the solid angle. The modules are read out at both ends by photomultipliers 
providing energy deposit and arrival time information. The readout segmen- 
tation provides the coordinates transverse to the fiber plane. The coordinate 
along the fibers is obtained by the difference between the arrival times of the 
signals at either end. Cells close in time and space are grouped into calorime- 
ter clusters. The energy and time resolutions are a E /E = 5.7%/ (GeV) 

and cr T = 54 ps/^E (GeV) © 50 ps, respectively. 

The KLOE trigger, described in Ref. 9, uses calorimeter and chamber in- 
formation. For this analysis, only the calorimeter signals are used. Two energy 
deposits above threshold (E > 50 MeV for the barrel and E > 150 MeV for 
the endcaps) are required. Recognition and rejection of cosmic-ray events is 
also performed at the trigger level. Events with two energy deposits above a 
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30 MeV threshold in two of the outermost calorimeter planes are rejected. 

During 2002 data taking, the maximum luminosity reached by DA$NE 
was 7.5 x 10 31 cm _2 s _1 , and in September 2002, DA$NE delivered 91.5 pb _1 . 
We collected data in 2001-2002 for an integrated luminosity L = 450 pb _1 . A 
total of 1.4 billion <fi mesons were produced, yielding 450 million Kg-K^ pairs. 
Assuming BR(K S — > 3tt°) = 1.9 x 10 -9 , ~1 signal event is expected to have 
been produced. 

The mean decay lengths of the K s and K L are \ s ~ 0.6 cm and \ L ~ 340 
cm at DA$NE. About 50% of K^s reach the calorimeter before decaying. 
The Kl interaction in the calorimeter ( u Kl crash" ) is identified by requiring 
a cluster with energy greater than 100 MeV that is not associated to any 
track and whose time corresponds to a velocity in the rest frame, (3*, of 
~ 0.2. The K L -crash provides a very clean K s tag. The average value of the 
center-of-mass energy, W, is obtained with a precision of 30 keV for each 100 
nb _1 running period (of duration ~ 1 hour) using large-angle Bhabha events. 
The value of W and the lf£,-crash cluster position allows us to establish, for 
each event, the trajectory of the Ks with an angular resolution of 1° and a 
momentum resolution better than 2 MeV. 

Because of its very short lifetime, the displacement of the K s from the 
decay position is negligible. We therefore identify as Ks decay photons neutral 
particles that travel with (3 = 1 from the interaction point (IP) to the EMC. 
Each cluster is required to satisfy the condition |T — R/c\ < min(3<7<r,2 ns), 
where T is the photon flight time and R the path length; ot also includes 
a contribution from the finite bunch length (2-3 cm), which introduces a 
dispersion in the collision time. 

In order to retain a large control sample for the background while preserv- 
ing high efficiency for the signal, we keep all photons satisfying E 1 > 7 MeV 
and | cos(0)| < 0.915. The photon detection efficiency is ~90% for E 1 = 20 
MeV, and reaches 100% above 70 MeV. The signal is searched for by requiring 
six prompt photons after tagging. 

The normalization is provided by counting the Ks — > 2tt° events in the 
same tagged sample. 

3 Monte Carlo simulation 

The response of the detector to the decay of interest and the various back- 
grounds is studied using the Monte Carlo (MC) program GEANFI [11]. GEAN- 
FI accounts for changes in machine operation and background conditions, 
following the machine conditions run by run, and has been calibrated with 
Bhabha scattering events and other processes. The response of the EMC to 
Kl interactions is not simulated but has been obtained from a large sample 
of ii'i-mesons tagged by identifying Ks — > 7r + 7r~ decays. This not only gives 
accurate representation of the EMC response to the K L crash, but also results 
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in an effective 40% increase in MC statistics. The fT^-crash efficiency cancels 
in the final 37r°/27r° ratio to better than 1% and we assign a 0.9% systematic 
error to the final result due to this source. 

Backgrounds are obtained from MC — > KgK^ events corresponding to an 
integrated luminosity L = 900 pb^ 1 . We also use a MC sample of <fi — > K + K~ 
events for L = 450 pb -1 and a MC sample of radiative decays for L = 2250 
pb _1 . A sample of ~ 340 000 K s — > 37r° MC events is used to obtain the signal 
efficiency. 



4 Photon counting for data and Monte Carlo 



To test how well the MC reproduces the observed photon multiplicity after 
tagging, we determine the fraction of events of given multiplicity, iV 7 = k, 
defined as F(k) = N ev (N^ = k)/ £? =3 N^N^ = i). As shown in Table 1, 
there is a significant discrepancy between data and Monte Carlo for events 
with multiplicity five and six. These samples are dominated by Kg — > 2tt° 





Data 2001 


MC 2001 


Data 2002 


MC 2002 


F(3) 


30.95±0.16 


30.31±0.11 


30.79±0.12 


30.06±0.08 


F(4) 


67.35±0.23 


67.93±0.17 


67.93±0.18 


68.15±0.12 


F(5) 


1.55±0.01 


1.80±0.01 


1.19±0.01 


1.66±0.01 


F(6) 


0.15±0.01 


0.14±0.01 


0.08±0.01 


0.13±0.01 



Table 1 

Measured values of F for data and Monte Carlo samples, in percent. 

decays plus additional clusters due either to shower fragmentation (split clus- 
ters) or the accidental coincidence of machine background photons (acciden- 
tal clusters). To understand this discrepancy, we have measured the proba- 
bility, Pa(1,2), of having one, or more than one, accidental cluster passing 
our selection by extrapolating the rates measured in an out-of-time window, 
(—68 < T < —14) ns, that is earlier than the bunch crossing. In Table 2, we 
list the average values of these probabilities. The observed discrepancy has 



F(K) 



Data 2001 MC 2001 Data 2002 MC 2002 



10 2 xP A (l) 0.75±0.30 1.03±0.16 0.38±0.17 0.89±0.08 

10 2 xP A (2) 0.14±0.05 0.16±0.03 0.07±0.02 0.10±0.03 

10 3 xP s (l) 3.6 ±0.2 3.8 ±0.3 3.7 ±0.2 3.3 ±0.1 

10 4 xP s (2) 1.5 ±0.4 1.5 ±0.3 0.9 ±0.2 1.7 ±0.2 

Table 2 

Measured values of the probabilities Pa and P$. 

been traced to an understood problem with the procedure for the selection of 
machine-background clusters. 

The MC-true fraction of events with a given multiplicity, /mc? is obtained 
by ignoring clusters due to machine background and counting at most one 
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cluster per simulated particle incident on the calorimeter. Using the fractions 
Jmc-i together with the values of Pa obtained as discussed above, we fit the 
observed F(k) distribution to get the probability for a cluster to generate frag- 
ments, P s (see Table 2). This fit accurately reproduces the observed fractions 
in the multiplicity bins five and six. More details on these measurements can 
be found in Ref. 12. The results of this study demonstrate the need for care- 
ful calibration of the background composition when comparing data and MC 
samples. 



5 Data analysis 

Kg — > 3ir° candidates consist of a Kl crash plus six photons. In our data 
sample of L = 450 pb -1 , we find 39 538 events, essentially all background. 
After removing background, we obtain the branching ratio by normalizing to 
the number of K$ — > 2ir° events. The latter are found by asking for three to 
five prompt photons plus the -fT^-crash. 

According to the MC, the six-photon sample is dominated (95%) by K s — > 
2tt° decays plus two additional photon clusters. These clusters are due to 
fragmented or split showers (2S, 1S+1A, 34%) and to accidental photons from 
machine background (2A, 64%). About 2% of the background events are due 
to false i^-crash tags from — » K$Kl — > 7r + 7r~, 3tt° events. In such events, 
charged pions from K s decays interact in the low-beta insertion quadrupoles 1 , 
ultimately simulating the f^-crash signal, while K L decays close to the IP 
produce six photons. Similarly, — > K + K~ events give a false signal (~ 1%), 
as well as <fi — > 777 — > 37r°7 events (~ 0.3%). The cuts described in the following 
make these contaminations negligible. 

To reduce the background, we first perform a kinematic fit with 11 con- 
straints: energy and momentum conservation, the kaon mass and the velocity 
of the six photons. The x 2 distribution of the fit to data and MC background 
is shown in Fig. 1. In the same plot, we also show the expected shape for signal 
events. Cutting at a reasonable x 2 value (x 2 /ll < 3) retains 71% of the signal 
while considerably reducing the background from false i^-crash events (33%), 
in which the direction of the Ks and Kl are not correlated. However, this cut 
is not as effective on the 2S, 2A background, due to the soft energy spectrum 
of the fake clusters. In order to gain rejection power over the background for 
events with split and accidental clusters, we look at the correlation between 
the following two x 2 -like estimators: 

• ( 2 , defined as 

_Amj Ami , (E Ks - E^ 7 J 2 (P? 3 -12iP?) 2 , 
< < °l ° 2 p 



The first quadrupoles are located approximately 45 cm on either side of the IP. 



6 




Fig. 1. Distribution of \ 2 f° r the tagged six-photon sample for data (points), MC 
background (solid line), and 10 000 events of MC signal (dashed line). 

(P^-H^f (Pf*-E,P?) 2 (7T-Q 2 



selecting the four out of six photons that provide the best kinematic agree- 
ment with the Kg — > 27r° decay hypothesis. This variable is quite insensitive 
to fake clusters. It is constructed using the two values of Am = m; — m n o 
(where m.j is the invariant mass of a photon pair), the opening angle between 
7r°'s in the K$ rest frame, and 4-momentum conservation. The resolutions 
on these quantities have been evaluated using a control sample of events 
with a f^-crash and four prompt photons. 
(s, defined as 



Am\ Am\ Am 2 

r 2 /t-2 /t-2 



3 



<T O" (T 

r 0v 



where the pairing of the six photons into n s is performed by minimizing 
this variable. (3 is close to zero for a ii's — > 3n° event and large for six- 
photon background events. 

For each estimator, the photon pairing with smallest £ value is kept. Figure 2a 
shows the distribution of events in the C2-C3 plane for the MC background. 
Most of the events are concentrated at low values of (2, as expected for K$ — > 
27T° events plus some additional isolated energy deposits in the EMC. A clear 
signal/background separation is achieved as can be seen by comparing the 
background and signal distributions in Figs. 2a and 2b. We subdivide the (2- 
(3 plane into the six regions Bl, B2, B3, B4, B5, and S as indicated in Fig. 
2a. Region S, with the largest signal-to-background value, is the "signal" box. 

The scatter plot in the C2-C3 plane for the data is shown in Fig. 3a. Our 
MC simulation does not accurately reproduce the absolute number of 2S and 
2A background events. This is also true of the predicted number of false Kl- 
crash events. However, the simulation does describe the kinematical properties 
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of these events quite well. The two-dimensional C2-C3 distribution allows us to 
calibrate the contributions from the different backgrounds. The MC shapes 
for each of the three categories are shown in Figs. 3b-d. We perform a binned 
likelihood fit of a linear combination of these shapes to the data, excluding 
the signal-box region. From the fit we find the composition of the six-photon 
sample to be (37.9±1.0)%, (57.4±1.3)%, and (4.7±0.3)% for the 2S, 2A, and 
false A^-crash categories, respectively. 

As a check, we compare data and MC for the projected distribution in 
£3 for the three bands in £ 2 , as shown in Figs. 4a-b. Excellent agreement is 
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Fig. 2. Scatter plot of (2 vs. £3 plane for the tagged six-photon sample: a) MC 
background, b) MC signal. 




Fig. 3. Scatter plots of C2 vs. £3 for the tagged six-photon sample: data (a), MC 
sample with two split clusters (b), two accidental clusters (c), and false i^L-crash 
events (d). 
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Bl 


B2 


S 


B3 


B5 


B4 


data 
MC 


452 ± 21 
419 ± 19 


10132 ± 101 
9978 ± 104 


282 ± 17 
282 ± 13 


5037 ± 71 
4816 ± 43 


326 ± 18 
380 ± 10 


22309 ± 149 
22682 ± 190 



Table 3 

Comparison between data and MC expectations in the different regions of the C2-C3 
plane for the entire sample with a -fT^-crash and six prompt photons. The boxes are 
defined as in Fig. 2a. 
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Fig. 4. Distributions in £3 for the tagged six-photon sample. Plots on the left are 
for events in the central band in £2; plots on the right are for events in all other 
regions of the plane. The plots in the top row are for the entire sample, before any 
cuts are made. The plots in the bottom row are after the application of the track 
veto. In all cases, black points represent data; solid line represents MC. 



observed. The large peak at low values of £3 in the central band is due to the 
false i\£-crash events. As a final test, we compare data and MC in the signal 
box and the five surrounding control regions. The agreement is better than 
10% in all regions, as seen from Table 3. 

Although cutting on x 2 substantially suppresses the false fr^-crash back- 
ground, we reduce this background to a negligible level by vetoing events with 
tracks coming from the IP. This effectively eliminates events in which the false 
Ifi-crash is due to a K s — > 7r + 7r~ decay with the pion secondaries interacting 
in the quadrupoles. The effect on the signal region can be appreciated by com- 
parison of Figs. 4a and 4c. Moreover, in order to improve the quality of the 
photon selection using ( 2 , we cut on the variable AE = (m^/2 — J2Ei)/cr E , 
where % = 1-4 stands for the four chosen photons in the (2 estimator and 
oe is the appropriate resolution. For K s — > 2n° decays plus two background 
clusters, we expect AE ~ 0, while for Ks — > 7r 7r°7r°, AE > ra^jo E . 
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25 50 75 100 50 100 



Fig. 5. Distributions of £3 for the central band 12. 1< £2 <60 (a), the side-bands 
C2 <12.1, C2 >60 (b), after all cuts. Points represent data, solid line MC. 





Bl 


B2 


S 




B3 


B4 


B5 


data 





4 ± 2 


2.0 ± 


1.4 


520 ± 23 


3 ± 2 


326 ± 18 


MC 





3.2 ± 0.8 


3.1 ± 


0.8 


447 ± 10 


2.5 ± 0.8 


389 ± 10 



Table 4 

Same as Table 3, after all cuts. 

Before opening the signal box, we refine our cuts on x 2 , ( 2 , C3; an d Ai±7 
using the optimization procedure described in Ref. 13. We end up choosing 
X 2 < 40.4 and AE >1.7. The signal box is defined by 12.1 < C2 < 60 and (3 < 
4.6. 

Figures 5a and 5b show the (3 distributions for the central band and the 
sidebands in ( 2 - 

In Table 4, we also list the number of events obtained in each of the six 
regions of the C2-C3 plane at this final stage of the analysis. In Figs. 6a-b 
we show the C2-C3 scatter plots for data and Monte Carlo. The rectangular 
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Fig. 6. Distribution of C2 vs (3 after cuts: MC background 900 pb _1 (a), data 450 
pb- 1 (b). 

region illustrates the boundaries of the optimized signal box. Seventeen MC 
events are counted in this region before applying the data-MC scale factors 
resulting from the calibration procedure described above. Contributions to 
the scale factors include the fact that the simulated integrated luminosity is 
greater than that for the data set (x2), the increased -K^-crash efficiency in 
the simulation (xl.4), and the increased probability of having accidental or 
split clusters in the simulation (on average, ~xl.9). 
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The selection efficiency at each step of the analysis has been studied using 
the MC. After tagging, the efficiency for the six-photon selection is (47.8 ± 
0.1 stat )%. Including all cuts, we estimate a total efficiency of e 37r = (24.4 ± 
0-l s tat)%- At the end of the analysis chain, we have two candidates with an 
expected background of B cxp = 3.13 ± 0.82 stat . 

In the same tagged sample, we also count events with photon multiplicities 
of three, four, or five. The corresponding efficiency is (91.8±0.2 sta t)% for Ks — > 
7r°7T events. The residual background contamination is estimated to be (0.77 
±0.24 stat+sys )% and (0.65±0.10 stat+sys )% in the 2001 and 2002 running periods 
respectively. Subtracting the background and correcting for the efficiency, we 
count 3.78 x 10 7 K s — > 7r°7r° events. We use this number to normalize the 
number of signal events when obtaining the branching ratio. 

6 Systematic uncertainties 

Systematics arise from uncertainties in estimation of the acceptance, back- 
grounds, and the analysis efficiency. The evaluation of the systematic uncer- 
tainties is described in detail in Ref. 14. 

Concerning the acceptance of the event selection for both the 2tc° and 37r° 
samples, we estimate the systematic errors in photon counting by comparing 
data and MC values for the P A and P s probabilities described above. The 
photon reconstruction efficiency for both data and MC is evaluated using a 
large sample of <fi — > 7r + 7r~7r°, 7r°^77events. The momentum of one of the 
photons is estimated from tracking information and position of the other clus- 
ter. The candidate photon is searched for within a search cone. The efficiency 
is parameterized as a function of the photon energy. Systematics related to 
this correction are obtained from the variation of the efficiency as a function 
of the width of the search cone. The results are listed in Table 5 under the 
heading cluster. The total uncertainty is smaller for the normalization sample 





Aa/a (Ks —> 


2tt°) 


Aa/a (Ks 




3vr°) 


Cluster 


0.16 % 




0.70 ' 


% 




Trigger 


0.08 % 




0.08 ' 


% 




Background filter 


0.20 % 




0.08 ' 


% 




Total 


0.27 % 




0.71 ' 


% 





Table 5 

Systematic acceptance uncertainties, Aa, for the 2ir° and 3ir° event selection crite- 
ria. 

since an inclusive selection criterion is used in this case. 

The normalization sample also suffers a small (0.4%) loss due to the use 
of a filter during data reconstruction to reject cosmic rays, Bhabha fragments 
from the low-beta quadrupole, and machine background events. This loss is 
estimated using the MC. We correct for it and add a 0.2% systematic error to 
the selection efficiency. The trigger and cosmic-ray veto efficiencies have been 
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estimated with data for the normalization sample and extrapolated by MC 
to the signal sample. These efficiencies are very close to unity and the related 
systematics are negligible. 

For the tagged six-photon sample, we have investigated the uncertainties 
related to the estimate of the background in the signal box after all cuts, B cxp . 
We have first considered the error related to the calibration of the MC back- 
ground composition by propagating the errors on the scale factors obtained 
from the fit. This corresponds to a relative error of 2.4% on B cxp . Moreover, 
we have investigated the extent to which the track-veto efficiency influences 
the residual false -K^-crash contamination. To do so, we examine the data-MC 
ratio, Ftp, of the sidebands in (3* for events rejected by this veto, since for true 
K^s (3* peaks at ~ 0.2 while false i^-crashes are broadly distributed in (3*. 
We obtain R p = 1.10 ± 0.01. Knowing that in the MC only 24% of the fakes 
survive the veto, we find a fractional error of 32% on the fake background. 
Since false i^-crash events account for 15% of the total background, the error 
on B cxp from data-MC differences in the track veto efficiency is 4.6%. 

A control sample of K s — > 2n° with four prompt photons has been used to 
compare the energy scale and resolution of the calorimeter in data and in the 
MC. The distributions of the m and AE variables have also been compared 
by fitting them with Gaussians. By varying the mass and energy resolution 
by ±lo" in the definitions of (2 and (3, we observe a relative change of 6.6% 
in the background estimate. Similarly, correcting for small differences in the 
energy scale for data and MC, we derive a systematic uncertainty of 6.7% on 
B exp . 

Finally, we have tested the effect of the cut on \ 2 by constructing the 
ratio between the cumulative distributions for data and MC. An error of 5% is 
obtained. A summary of all the systematic errors on the background estimate is 
given in Table 6. Adding in quadrature all sources we obtain a total systematic 





AB cxp /B cxp 




Background composition 


2.4% 




Track veto 


4.8% 


0.2% 


Energy resolution 


6.6% 


0.5% 


Energy scale 


6.7% 


1.0% 


x 2 


5.0% 


1.8% 


Total 


11.5% 


2.1% 



Table 6 

Systematic uncertainties on the expected background and analysis efficiency, e ana . 

error of 12% on the background estimate. 

To determine the systematics related to the analysis cuts for the signal, we 
have first evaluated the effect of the track veto. Using the MC signal sample, 
we estimate a vetoed event fraction of (3.7 ± 0.1)%. The data-MC ratio of the 
cumulative distributions for the track-vetoed events in the tagged six-photon 
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sample is Rtv = 1.06 ±0.03, which translates into a 0.2% systematic error on 
the track-veto efficiency. 

Because of the similarity of the \ 2 distributions for the tagged four- and 
six-photon samples, as confirmed by MC studies, an estimate of the systematic 
error associated with the application of the x 2 cu ^ can be obtained from the 
data-MC comparison of the cumulative x 2 distributions for the four-photon 
sample. The systematic error arising from data-MC discrepancies in the \ 2 
distribution is estimated to be 1.8% by this comparison. 

Moreover, the efficiency changes related to differences between the calorime- 
ter resolution and energy scale for data and MC events have been studied in a 
manner similar to that previously described for the evaluation of the systemat- 
ics on the background. The systematic uncertainties on the analysis efficiency 
are summarized in Table 6. Adding all sources in quadrature we quote a total 
systematic error of 2.1% on the estimate of the analysis efficiency. 

7 Results 

At the end of the analysis, we find 2 events in the signal box with an estimated 
background of B exp = 3.13 ± 0.82 sta t ± 0.37 syst . To derive an upper limit on 
the number of signal counts, we build the background probability distribution 
function, taking into account our finite MC statistics and the uncertainties on 
the MC calibration factors. This function is folded with a Gaussian of width 
equivalent to the entire systematic uncertainty on the background. Using the 
Neyman construction described in Ref. 15, we limit the number of Kg — > 37r° 
decays observed to 3.45 at 90% C.L., with a total reconstruction efficiency of 
(24.36 ± 0.11 s t a t ± 0.57 sys )%. In the same tagged sample, we count 3.78 x 10 7 
Kg — > 7r°7T events. This number is used for normalization. Finally, using the 
value BR(K S -> tt°tt°) = 0.3105 ± 0.0014 [16] we obtain: 

BR(K S -> 3tt°) < 1.2 x 10~ 7 at 90% C.L., (1) 



which represents an improvement by a factor of ~6 with respect to the best 
previous limit [4], and by a factor of 100 with respect to the best limit obtained 
with a direct search [2]. 

The limit on the BR can be directly translated into a limit on |r?oool : 



l^ooo I = 



A(K S - 


+ 3tt°) 






A(K L - 


^3tt°) 










r L BR(K S - 


^3tt°) 




= \ 


r s BR(K L - 


-> 3tt°) 



(2) 



< 0.018 at 90% C.L. 



This result describes a circle of radius 0.018 centered at zero in the J^ooo); 
Q^ooo) plane and represents a limit 2.5 times smaller than the result of Ref. 4. 
As follows from the discussion in that reference, our result confirms that the 
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sensitivity of the CPT test from unitarity is now limited by the uncertainty 
on r] + _. 
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